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ABSTRACT 

We investigate the behavior of the magnetic Prandtl number (ratio of mi- 
croscopic viscosity to resistivity) for accretion sources. Generally this number is 
very small in standard accretion disk models, but can become larger than unity 
within ~ 50 Schwarzschild radii of the central mass. Recent numerical investiga- 
tions suggest a marked dependence of the level of MHD turbulence on the value 
of the Prandtl number. Hence, black hole and neutron star accretors, i.e. com- 
pact X-ray sources, are affected. The astrophysical consequences of this could 
be significant, including a possible route to understanding the mysterious state 
changes that have long characterized these sources. 

Subject headings: accretion, accretion disks; black hole physics; magnetic fields; 
MHD; turbulence 



1. Introduction 

Magnetohydrodynamic (MHD) turbulence differs from ordinary hydrodynamic turbu- 
lence in at least one very important respect: whereas the latter generally has only one 
dissipative scale (viscous), MHD turbulence has two (viscous and resistive). This raises the 
question of whether the classical Kolmogorov picture, in which the large scale energetics 
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of the turbulenct cascade is essentially independent of small scale dissipation physics, re- 
mains valid in MHD turbulence. And if it is not valid, what are some possible astrophysical 
consequences? 

The dimensionless ratio of the kinematic viscosity v to the electrical resistivity r\ is 
known as the magnetic Prandtl number, Pm. Balbus & Hawley (1998) suggested that even 
if both the viscous and resistive dissipation scales are very small, the saturation level of the 
MHD turbulence produced by the magnetorotational instability (MRI) should be sensitive to 
Pm, at least in the neighborhood of Pm ~ 1. Their argument was as follows. If Pm <C 1, the 
resistive scale is much larger than the viscous scale. Assuming that the velocity fluctuations 
do not greatly exceed the Alfvenic fluctuations at the resistive scale, viscous stresses on 
the resistive scale would be small. This would mean that even if relatively large velocity 
gradients accompanied magnetic dissipation, these gradients would not produce stresses that 
would interfere with the dynamics of the field reconnection. Since Lorentz forces drive the 
MRI, the dissipation of the magnetic field is an important regulatory mechanism for the 
saturation level of the turbulence. On the other hand, if Pm ^> 1, and the viscous scale 
were significantly larger than the resistive, then the resulting dynamical stresses would likely 
be relatively large when the magnetic field is dissipated at the small resistive scale. (This 
assumes that significant velocity fluctuations accompany small scale reconnection. Such 
fluctuations would be heavily damped at the resistive scale, and any reconnection would 
have to be very slow.) These stresses would then interfere with the field reconnection and 
dissipation, leading to a build-up of magnetic energy that cascades upwards, back to larger 
scales (Brandenburg 2001). 

Balbus & Hawley (1998) were motivated by the possibility that the turbulent properties 
of accretion disks might be different in the regimes Pm ^> 1 and Pm < 1. At the time, 
direct numerical simulation of flows with different Prandtl numbers was very difficult, and 
these authors attempted only the crudest of tests by varying the level of artificial viscosity 
in the ZEUS MHD code at fixed resolution. These preliminary experiments did, however, 
show a higher level of saturation for a larger viscosity. Since this is an example in which 
increasing a dissipation coefficient actually raised the level of turbulent activity, it was a 
noteworthy result. 

A decade on, it is possible to do much better. There is a definite sensitivity to Pm 
in numerical simulations of MHD turbulence. A Pm dependence has been observed for a 
number of years now in stirred magnetic turbulence (Schekochihin et al. 2004, 2005). For 
example, fluctuation dynamos at large Pm were found numerically by Schekochihin et al. 
(2004), but until recently there was some question as to whether a low Pm fluctuation dynamo 
even existed; this has now been answered affirmatively (Iskakov et al. 2007; Schekochihin 



- 3- 



et al. 2007). For astrophysical accretion flows, MRI calculations are of direct interest, and 
the last year has seen the first Pm studies in "shearing box" simulations. Zero mean field 
calculations have been carried out by Fromang et al. (2007), while Lesur & Longaretti (2007) 
studied a mean vertical field. In the latter investigation, the radial-azimuthal component 
of the stress tensor behaved linearly over the range 0.12 < Pm < 8 with no apparent 
sign of approaching an asymptote (Iskakov et al. [2007] did, however, appear to be reaching 
saturation levels in some of their driven turbulence runs) . With extensive numerical evidence 
of a Pm dependence in MHD turbulence, a natural question to raise is what is the behavior 
of Pm in classical accretion disk models? In particular, is there a transition from Pm <C 1 
to Pm 3> 1 in phenomenological models that have been used to model AGN and compact 
X-ray sources? 

In this paper, we examine the magnetic Prandtl number behavior of classical a models. 
In fact, the only feature of these models that is important for our purposes is that the free 
energy of differential rotation be locally dissipated — a variable a parameter, for example, 
would hardly change our conclusions at all. Throughout the regime of interest, the disk 
is fully ionized and collision dominated (see §2 below), so that the Spitzer (1962) values 
for the resistivity and viscosity are appropriate. Our principal finding is that generally 
Pm <C 1 nearly everywhere in classical a models, with one robust and important exception: 
on scales less than ~ 100 Schwarzschild radii in black hole and neutron star disks. It is 
extremely tempting to associate this Prandtl number transition with a physical transition 
in the properties of the accretion flow, here motivated by "first principle" physics. Further 
discussion of this point is presented below. 

An outline of the paper is as follows. §2 presents preliminary estimates of important 
parameter regimes. §3 is the heart of the paper, in which we calculate the behavior of Pm 
in classical a disk models. Transitions from low Pm to high Pm regions occur only in disks 
around black holes and neutron stars. Finally, §4 is a discussion of the possible astrophysical 
consequences of having both high Pm and low Pm regions in the same disk. It is argued 
that high and low X-ray states (e.g. McClintock & Remillard 2006) may be related to an 
unstable interface between Pm < 1 and Pm > 1 regions of the disk. 

2. Preliminaries 

The magnetic Prandtl number is not a standard parameter of accretion theory, so let us 
begin with a brief orientation in the temperature-density parameter space. Throughout this 
work, the fiducial disk plasma is taken to be a mixture of 90% hydrogen and 10% helium (by 
number). Following the discussion in Spitzer (pp. 138-9), we estimate an averaged resistivity 
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of such a fully ionized gas as 

5.55 x 10 n lnA eH 2 -i 

V = f3f2 Cm S > 

where T is the temperature in Kelvins, and A c h is the Coulomb logarithm for electron-proton 
scattering. (Modifications in the logarithm due to electron-helium scattering, here a minor 
effect, are ignored.) 

The kinematic viscosity of the same gas is estimated to be 

1.6 x 1(T 15 T 5 / 2 o -i 

v = — cm 2 s" 1 , 2 

plnA HH 

where p is the mass density and lnA H H is the Coulomb logarithm for scattering of protons 
by protons. (See Appendix for a derivation of these results and a discussion of the Coulomb 
logarithms.) This gives a magnetic Prandtl number of 

Pm = 2.9 x KT 27 T , . (3) 

plnA eH m A H h 

The two logarithms differ from one another for temperatures in excess of 4.2 x 10 5 K (see 
Appendix). If / is the product of the two Coulomb logarithms normalized to a nominal value 
of 40, 

T \ 4 /l0 22 cm- 3 \ / T \ 4 /l0 14 cm- 3. 



Pm V4.2x10 6 kJ \ ln H ) V4-2x10 4 kJ \ I n H J' ^ 

where n# is the number density of hydrogen atoms. The two last forms that are given 
for Pm are convenient for applications to a 10 M (binary) and 10 8 M Q (AGN) black hole, 
respectively. 

Finally, it is required to justify quantitatively the statement in the Introduction that 
the disk plasma is collisional near the transtion point Pm = 1. We shall refer to a plasma 
as "dilute" (as opposed to collisional) if the product of the ion cyclotron frequency u ci and 
the ion-ion collision time t ci is greater than one. The proton cyclotron frequency may be 
written 

-'Off 

where we have introduced the plasma (5 parameter, the ratio of the gas to magnetic pressure. 
For a gas of cosmic abundances, 



w„ = 8.6 x lO- PiH > ( 5 ) 
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In a dilute plasma, it is not appropriate to use the Spitzer (1962) form of the viscosity, as 
we have done above, hence we need a numerical estimate of w C it C i. (It is also not strictly 
correct to use the Spitzer resistivity, but the correction here is relatively minor.) Following 
the prescription set forth in the Appendix (divide the nominal proton-proton collision time 
by a factor of 1.5 to include the effects of proton-helium collisons), we obtain 

(cu ct t ci ) 2 ~ 9 x 1(T 7 -^-. (7) 

n H fj 

where we have taken the relevant Coulomb logarithm to be 7. This should be compared 
directly to Pm ~3x 10~ 5 T 4 /n# from equation ([3]). For a given value of (3, the temperature 
and density dependence of co^i and Pm are the same. What is more, we are concerned 
in this work with weakly magnetized plasmas, (3 > 1, and generally ,8 > 1. Therefore, at 
the threshold Pm = 1, the plasma is never dilute, and the collisional regime is valid. Note, 
however, that once into the large Pm regime, substantial heating and magnetic field growth 
may lead to a dilute plasma phase, or perhaps even to a fully collisionless phase in which 
the fluid approximation itself breaks down. 



3. Analysis 

Our goal is a simple one: we wish to follow the behavior of Pm with disk radius in a 
standard a model, in effect testing such models for self-consistency. If most of the energy 
extracted from differential rotation is locally dissipated, the basic a scalings are probably 
robust. This is particularly true if the problem is framed to minimize any possible explicit 
dependence upon a of the temperature and density, as we have done. Then, even if in real 
disks it is not a very good approximation to treat a as a constant, its variability is not 
crucially important for the scaling laws. 



3.1. Pm behavior in a models 

Our starting point is the Kramers opacity disk model of Frank, King, and Raine (2002). 
The density in the midplane is 

p = 3.1 x lO-^-^M^^M/Mo) 5 / 8 ^- 15 / 8 ? 11 / 20 gmcm" 3 , (8) 

where M\§ is the mass accretion rate in units of 10 16 g s _1 , M/M & is the central mass in 
solar units, Rio is the cylindrical radius R in units of 10 cm., and q = 1 - (R^R) 1 ' 2 . The 
quantity R* is a fiducial radius at which the stress is taken to vanish (the "inner edge"), 
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but in practice we shall assume that R ^> R*, and hence that q is unity. The midplane 
temperature is given by Frank et al. (2002) as 

T = 1.4 x lO 4 a- 1 / 5 Aff/ lo (M/M ) 1 / 4 J R lo - 3 / 4 g 3/io K (9) 
This leads to a Prandtl number of 

Pm = 9.0 x 10- 5 l-^ 10 MS /x (M/M Q )^R 10 - 9 V^. (10) 

Typical disk Prandtl numbers are therefore very small, and insensitive to scaling with a. 
Transitions from low to high Pm, if they occur at all, will occur in the inner disk regions. 

Let us calculate R cr , the critical radius at which Pm = I. Here, it will suffice to set 
q = 1 (R 3> it!*); a more accurate numerical calculation (described below) certainly justifies 
this. With Pm = 1, we find 

R cr = 2.5 x 10 6 r 8 / 9 a-V 45 M™ /45 (M/M Q y/ 3 cm. (11) 

The region of interest is evidently on scales of tens of Schwarzschild radii (Rs)- With 
R s = 2GM/c 2 , this becomes 

^ = ^a^Mlt^iM/M^H-^. (12) 

Rs 

Our final step is to scale the mass accretion rate with M. If we assume that the source 
luminosity L is a fraction e of Mc 2 and a fraction S of the Eddington luminosity 

L Edd = 1.26 x 1O 38 (M/M ) erg s" 1 , 

then 

^ = 59 (a_ 2 M/M )" 4/45 (5/e) 26/45 r 8/9 (13) 
Rs 

where a_ 2 is a in units of 0.01. The ratio S/e is just the mass accretion rate measured in units 
of the Eddington value Mem = L E dd/c 2 . This shows that the critical radius at which the 
Prandtl number transition occurs, when measured in units of Rs, is remarkably insensitive 
to the central mass. In general we find that R cr varies roughly between 10 and 100 Rs- 
In principle, the low Pm region could in some cases extend all the way down to 2 — 3Rs, 
particularly for larger AGN masses. Iron line observations of, for example, the well-studied 
Seyfert galaxy MCG-60-30-15, (Fabian et al. 2002) suggest the presence of an ordinary 
Keplerian-like disk down to 3Rs, and the Pm transition hypothesis must accommodate this: 
no transition should also be a possibility. 
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3.2. Pm behavior in numerical a models 

The result of the previous section neglects radiation pressure and electron scattering 
contributions to the opacity. In particular, the radiation to gas pressure ratio is easily 
calculated. With q — 1, 

^ = 5x lO-V' 10 ^ 2 °(M/M e ) 1 / 8 R 1Q - 3/8 . (14) 
(This differs from equation [5.56] in Frank et al. [2002].) At R = R cr , 

2/15 



(15) 



This varies between a 10% and an order unity effect for applications of interest. To ensure 
that radiative corrections do not alter the basic conclusion of the existence of a crtical Pm 
transition radius under nominal conditions, we have adapted the disk code of Terquem & 
Papaloizou (1999) to construct more detailed a models. Both radiation pressure and electron 
scattering opacity were included. We find that the essential qualitative features of equation 
(TIB"]) remain intact, though radiative effects do alter the scalings somewhat. We focus on two 
central masses, one a source of 10M Q (representative of an X-ray binary), the other 1O 8 M , 
which is representative of an AGN. The Prandtl number behavior for each of these cases for 
several different values of a, but at a fixed accretion rate (0.1M Edd ), is shown in figure (1). 
The two cases are very similar. Starting with a standard Keplerian a disk, these black hole 
accretion sources seem to make a transition from low Pm to high Pm at a typical value of 
~ 50^5- 

Figures (2) and (3) show Pm plots as meridional slices. A central mass of 10M Q is 
assumed for figure (2), while figure (3) corresponds to 10 8 M Q . In each figure, the left and 
right diagrams correspond respectively to M / M Edd = 0.001, 1. We have used a = 0.05. At 
higher accretion rates, the Pm > 1 region can be extensive; on the other hand, if M / Mem is 
sufficiently small, the flow can have Pm < 1 down to the marginally stable orbit R = 3Rs- 



4. Discussion 

The findings of the previous section show that if a models are even qualitatively correct 
in their scalings, only black holes and neutron star accretion disks, i.e., classical X-ray 
sources, will have regions with Pm < 1 and Pm > 1. If, as we would argue, there is 
a physical difference in the saturated state of MHD turbulence in these two regimes, it 
should be reflected in the astrophysical behavior manifested by X-ray sources. We tentatively 
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Fig. 1. — Behavior of Pm for 1OM (left) and 10 8 M Q (right) black holes for several different a 
values as a function of disk radius Rj Rs, where Rs is the Schwarzschild radius. The accretion 
rate is taken to be OAM^dd- Calculations were carried out using the code of Terquem & 
Papaloizou (1999) including radiation pressure and electron scattering corrections to the 
opacity. Note the insensitivity of the results to both a and the central mass. 




-100 -50 50 100 -600 -400 -200 200 400 600 

R/Rs R/Rs 



1e-3 1e-2 le-1 I.eO 1e1 1e2 1 e3 1e-3 1e-2 ls-1 I.eO 1e1 1e2 1 e3 

Magnetic Prandtl Magnetic prandt 

Fig. 2. — Regions of Pm for a 10M Q black hole accretion disk. Left diagram corresponds to 
M = O.OlMEdd, right to M = MEdd- In both cases, a = 0.05. At large accretion rates, the 
high Pm regions can be quite extended. (Note change of radial scale.) 
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suggest that the principal hard high states and low soft states associated with these sources is 
related respectively to the relative radiative dominance of the Pm > 1 and Pm < 1 accretion 
regions. In this discussion we will outline arguments that are suggestive, but as yet far from 
conclusive, of this. They are meant to spur further numerical investigation in what could 
prove to be an interesting direction. 

The results of several independent numerical simulations of MHD turbulence, both 
forced and shear-driven, appear to indicate that if Pm > 1, field dissipation becomes more 
inefficient, apparently because viscous stresses make the resistive scale less accessible (Fro- 
mang et al. 2007, Lesur & Longaretti 2007, Iskakov et al. 2007). If field dissipation is 
inefficient, the most likely scenario is that the field will initially build up on the viscous 



scale, but ultimately cascade upward to larger scales ([Brandenburg II2001I ). In a disk, the 
growing magnetic field would drive the MRI more vigorously until ultimately — and "ulti- 
mately" may in fact be rather rapid — the field is of order thermal. At this stage further 
MRI development is likely to be halted. 

The effective absence of resistivity of course does not mean that dissipation is absent; 
dissipative heating will still be present in the form of viscous heating. Note that the domi- 
nance of the resistive scale in Pm < 1 turbulence means that the electrons are directly heated 
(assuming that classical Spitzer resistivity applies), whereas the ions are directly heated in 
viscosity dominated Pm > 1 turbulence. The need for the dominance of ion heating in low 
luminosity black hole accretion is by no means a new idea (e.g. Narayan & Yi 1995), but 
placing it within the Prandtl number framework lends mutual support to this current work 
and to what has become the standard picture. In addition, the heating of a Pm > 1 mag- 
netized plasma may be very vigorous — unlike ohmic resistivity, viscous thermalization does 
not destroy the current sources. 

Conditions for a thermal runaway are present: at constant pressure, Pm oc T 5 . Since 
Pm is an increasing function of temperature, a little heating would tip Pm ~ 1 accretion 
towards the direction of Pm ^> 1 accretion. This would mean yet greater heating, following 
the numerical lead that large Pm turbulence is characterized by greater fluctuation levels. 
But this argument works in both directions, cooling as well as heating. A formal boundary 
between high and low Prandlt number regions would, in this view, be unstable because of 
the dependence of Pm upon T. This can be investigated by direct simulation. We speculate 
that the Pm < 1 region is a standard a disk and responsible for soft thermal emission; the 
Pm > 1 region corresponds to lower density, hotter accretion. Although such a component 
has been regarded as essential for understanding the X-ray spectra of black hole accretion 
sources, the reason that a disk would suddenly make the transition from one type of flow to 
the other has always been unclear. Prandtl number considerations may make this transition 
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less mysterious. 

The unstable boundary between high and low Prandtl number flow marks a fundamental 
change in the accretion flow, leading to a distinct inner accretion zone that dominates the 
hard tail of the X-ray spectrum. A better understanding of the interface will help to establish 
whether it is involved with transitions from one state to another. In a subsequent publication, 
we will present a technique to make this problem tractable and predictive. 

At this point the major gap in our scenario is the difference between the modest but 
rigorous numerical findings of a correlation between Pm with the amplitude of the turbulent 
stress, and the full blown thermal runaway that we envisage. That gap can begin to be filled 
with well-crafted numerical investigations of temperature-dependent dissipation coefficients 
in MRI turbulence. Little has yet been done along these lines, and it promises to be extremely 
challenging, particularly if the ideas presented here are qualitatively correct and a dilute or 
collisionless plasma appears. But preliminary investigations have already begun. 

We end by noting that in the paper introducing the MRI to the astrophysical community 
(Balbus & Hawley 1991), two possible nonlinear resolutions were envisioned. In one the 
field was limited to subthermal strengths by vigorous dissipation; in the other it grew to 
thermal levels and became "stiff." Subsequent numerical simulations seem to support the 
first outcome, but this may well have been because the codes used were not in the large Pm 
regime. Both scenarios might in fact be equally viable, the choice of direction being made 
by the Prandtl number of the turbulence. 
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Appendix: Collision time and viscosity estimates. 

Spitzer (1962) defines a "deflection time" t D for a test particle (denoted by subscript t) 
of mass m t , charge Z t (in units of e), and velocity w t to be scattered by Coulomb interactions 
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by a population of field particles (denoted by subscript /) of number density rif. It is given 
by 



2 3 



t D (t - /) - 8nnfe 4 Z 2 Z 2 F{lfWt) lnAt/ - ( 16 ) 



Here, 



(17) 

with m/ and 7/ refering respectively to the mass and temperature of the field particles. The 
function F(x) is 

F(*)=(l-^)erf W + ^ (18) 
where erf(x) denotes the standard error function 

erf(x) = 4= / e~ s2 ds. (19) 
V ff Jo 

The argument of the logarithm is 

1.5 /A; 3 T 2 T e \ 1/2 , , 

A t/ = —— l — 20 

Z t Z f e 6 \ irn e J 

where T t and T e refer to the test particle and electron temperature, respecitvely, and n e is 
the electron density. When the test particles are electrons, then for temperatures in excess 
of 4.2 x 10 5 K, an additional factor of (4.2 x 10 5 /T e ) 1/2 appears in the expression for A (a 
correction for quantum diffraction). 

In what follows, we shall always take a single temperature (T) fluid, and set w t equal 
to the rms test particle velocity, i.e., m t wl = 3kT. Then, 

t (t yf)= mTmf 12 m) 

87in f e 4 ZfZj F{^3m f /2m t ) In A tf ' 

and, with % denoting hydrogen number density, 

1.5 /A; 3 T 3 \ 1/2 _ 1.239 x 10 4 T 3 / 2 _ 1.131 x 10 4 T 3 / 2 

tf " W ) " z t z f ^ - z t z f ^ (22) 

with the additional diffraction correction of a factor of (4.2 x 10 5 /T) 1 / 2 needed for the case of 
electron test particles as noted above. For our cosmic gas, n e = 1.2n# under the assumption 
of fully ionized helium. If a single temperature prevails, then A t f = Aj t ; note that the time 
to(t — > f) does not have a similar symmetry between t and /. 
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As discussed in the text, representative values for T and % near the Prandtl number 
transition are T = 6 x 10 6 K and n H = 2 x 10 22 cm -3 . For these values, 

lnA HH = 7.07, lnA H He = 6.38, (23) 

and 

tD(H _ H) = L614 ( _™_) s . (24) 

The dynamical ion viscosity ^ of a fully ionized plasma is (Spitzer 1962): 

0.406m 1 / 2 (A;T) 5 / 2 , , 

where both the test and field particles are identified with ions of mass m and charge Z. 
Dimensionally, this takes the form 

Vv = C Vv pw 2 t D (26) 

where C Vv is a numerical constant (nominally but universally 1/3 in elementary modeling), 
p is the ion density, w 2 is the mean squared ion thermal velocity (3kT divided by the ion 
mass), and t D is the ion-ion deflection time. In considering a cosmic mixture of a 10% 
helium abundance fraction, one must take into account modifications to t D due to scattering 
of protons by He nuclei, in addition to the contribution to the viscous stress carried by these 
nuclei. Because of the sensitive dependence on atomic number Z, a relatively small amount 
of He could in principle make a significant contribution to rjy. Indeed, fully ionized metals 
at the level of a few per cent also make a contribution because of the Z scaling, but we shall 
ignore this here. Assuming that C Vv is the same for all species, an estimate for the cosmic 
abundance viscosity is then 

??(cosmic) = C ri (p H w 2 H t D (H) + p He w 2 He t D (He)) (27) 

where the deflection times are now given by 

1 - 1 + * „,= * ^fi + j#^y (28) 



t D (H) t D (H -> H) t D (H -> He) t D (H -> H) \ t D (H -> He) 

1 _ 1 1 _ 1 / t D (He -+H)\ 

t D (He) ~ t D (He -> H) + t D (He -> He) ~ t D (He -> H) \ + t D (He -»• He)) ' ^ 



Now, 



t D (H H) _ fn He \ x 4 J FWQ\ x /ta/W . 



t D (H^He) \n H J \F(y/TI) J V 1iiA hh 
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and 



t D (He H) _ (n H ,\ F(VU)\ ^ /In/W 



t D (He^He) \n H J lF(v^375)/ V lnA 



»-HcH 



In each of the above, the ratio of the Coulomb logarithms is about 0.9 across a wide range 
of densities and temperatures. Adopting this value, we find 

1 1.46 1 1.6 



t D (H) to(H — > H) ' t D (He) to(He — * H) 

In other words, the effects of test particles interacting with the 10% admixture of He results 
in roughly a 50% increase in the effective collision rate. At the level of accuracy with which 
we are concerned, we shall a deflection time shortening factor of 2/3 in both cases. The 
effective viscosity is then 



2 2 
-C v [p H w 2 H t D (H -> H) + p H eW 2 He t D {He -> if)] = -pz/f/ 



1 + ^£ x ^ x tz? ^ 6 ~" ^ 



p^ < t D (H^H) 

(33) 

where pi/# is the dynamical viscosity in a gas of pure hydrogen [y being the corresponding 
kinematic viscosity). The final deflection time ratio is 

w^^m = 2 x (V4) x [F(VU)/F(y/375)} x [In A HH / In A Hc h] ^ 0.925 (34) 

The final estimate for the cosmic abundance viscosity is 

2 _ T 5/2 

pi/(cosmic) = - x 1.09 x pv H = 1.6 x 10 15 - — gmcm x s 1 (35) 

3 hiA H H 
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Fig. 3. — Regions of Pm for a 1O 8 M black hole accretion disk. Left diagram corresponds to 
M = 0.01M E dd, right to M — M Edd . In both cases, a = 0.05. The high Pm regime is very 
limited in the low M case on the left. (Note change of radial scale.) 



